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ABSTRACT 

The problem of formation of the Ruderman-Sutherland type inner vacuum 
gap in neutron stars with ft ■ B < is considered. It is argued by means of the 
condition Ti/T s > 1 (where Tj is the critical temperature above which ^Fe ions 
will not be bound at the surface and T s is the actual temperature of the polar cap 
surface heated by the back-flow of relativistic electrons) that the inner vacuum 
gap can form, provided that the actual surface magnetic field is extremaly strong 
(B s > 10 13 G) and curved (TZ < 10 6 cm), irrespective of the value of dipolar 
component measured from the pulsar spin down rate. Calculations are carried 
out for pulsars with drifting subpulses and/or periodic intensity modulations, in 
which the existence of the quasi steady vacuum gap discharging via E x B drifting 
sparks is almost unavoidable. Using different pair-production mechanisms and 
different estimates of the cohesive energies of surface iron ions, we show that it 
is easier to form the vacuum gap controlled by the resonant inverse Compton 
scaterring seed photons than by the curvature radiation seed photons. 

1. Introduction 

The consecutive subpulses in a sequence of single pulses of a number of pulsars change 
phase systematically between adjacent pulses, forming apparent driftbands of the duration 
from several to a few tenths of pulse periods. The subpulse intensity is also systematically 
modulated along driftbands, typically increasing towards the pulse centre. In some pulsars, 
which can be associated with the central cut of the line-of-sight trajectory, only the periodic 
intensity modulation is observed, without any systematic change of subpulse phase. On 
the other hand, the clear subpulse driftbands are found in pulsars associated with grazing 
line-of-sight trajectories. These characteristics strongly suggest an interpretation of this 
phenomenon as a system of subpulse-associated beams rotating slowly around the magnetic 
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axis. Ruderman & Sutherland (1975, hereafter RS75) proposed a natural explanation of the 
origin of subpulse drift, which involved a number of isolated E x B drifting sparks discharging 
the quasi steady vacuum gap formed above the polar cap of the pulsar with ft ■ B < 0, in 
which ^Fe ions were strongly bound at the surface. Although the original idea of RS75 
associating the rotating sub-beams with the circulating sparks is still regarded as the best 
model of drifting subpulse phenomenon, their vacuum gap was later demonstrated to suffer 
from the so-called binding energy problem (for review see Abrahams & Shapiro 1991; Usov 
& Melrose 1995; Xu et al. 1999). In fact, the cohesive energies of ^Fe ions used by RS75 
proved largely overestimated and the inner vacuum gap envisioned by RS75 was impossible 
to form. However, it is worth emphasizing that RS75 considered the canonical surface 
dipolar magnetic fields with values determined from the pulsar spindown rate, although 
they implicitly assumed small radii of curvature 1Z ~ 10 6 cm required by the pair creation 
conditions, which is inconsistent with a purely dipolar field. Recently Gil & Mitra (2001; 
hereafter GM01) revisited the binding energy problem an found that the formation of a 
vacuum gap is in principle possible, although it requires an extremely strong non-dipolar 
surface magnetic field B s = bB d , where the coefficient & > 1 in a typical pulsar, B d = 
6.4 x 10 19 (PP) a5 G = 2 x 10 12 (PP_ 15 ) a5 G is the dipolar field at the pole, P is the pulsar 
period in seconds, P is the period derivative and P_i 5 = P/10~ 15 . 

In a superstrong surface magnetic field B s > 0.1B q , where B q = 4.414 x 10 13 G, the 
asymptotic approximation of Erber (1966) used by RS75 in derivation of the height of quasi 
steady vacuum gap is no longer valid. In fact, in such strong field the high energy Ef — hcu 
photons produce electron-positron pairs at or near the kinematic threshold hcu = 2mc 2 / sin 9, 
where sin# = h/TZ, h is the gap height, and 1Z = "T^IO 6 cm is the radius of curvature of 
surface magnetic field lines (e.g. Daugherty & Harding 1983), H is the Planck constant, c is 
the speed of light, m and e are the electron mass and charge, respectively. The vacuum gap 
formed under such conditions was called the Near Threshold Vacuum Gap (hereafter NTVG) 
by GM01. They considered two kinds of high energy seed photons dominating the e~e + pair 
production: the Curvature Radiation (CR) photons with energy fuv = (3/2)h'y 3 c/1Z (RS75; 
Zhang & Qiao 1996), and resonant Inverse Compton Scattering (ICS) photons with energy 
huj = 2'yheB s /mc (Zhang & Qiao 1996; Zhang et al. 1997), where 7 is a typical Lorentz factor 
of particles within the gap. The corresponding vacuum gap is called the Curvature Radia- 
tion dominated (CR-NTVG) and the Inverse Compton Scattering dominated (ICS-NTVG), 
respectively. GM01 estimated the characteristic heights of both CR-NTVG and ICS-NTVG. 
In this paper we further refine these estimates by including the general relativistic (GR) 
effects of inertial frame dragging (IFD) and considering the heat flow conditions within the 
thin uppermost surface layer of the polar cap. Moreover, we use a broader range of cohesive 
energies of surface ^Fe ions. The obtained VG models are applied to pulsars with drifting 
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subpulses and/or periodic intensity modulations, in which the presence of E x B drifting 
spark discharges seems almost unavoidable (Deshpande & Rankin 1999, 2001; Vivekanand 
& Joshi 1999, and Gil & Sendyk 2000). 



2. Near threshold vacuum gap formation 

If the cohesive energy of ^Fe ions is large enough to prevent them from thermionic (this 
section) or field emission (Appendix A), a vacuum gap forms right above the polar cap with 
a characteristic radius r p = 5-°- 5 io 4 P -0 - 5 cm (e.g. Gil & Sendyk 2000). Jessner et al. (2001) 
pointed out that GR effect of IFD (Muslimov & Tsygan 1992; Muslimov & Harding 1997) 
should affect the RS75 type models with a vacuum gap above the polar cap. Although Jessner 
et al. (2001) did not investigate the problem, they implicitly suggested that the electric fields 
distorted by GR effects make formation of "starved" inner magnetospheric regions even more 
difficult than in the flat space case. However, Zhang et al. (2000) demonstrated that although 
GR-IFD effect is small, it nevertheless slightly helps formation of VG above the polar caps. 
In other words, GR modified potential drop within VG is slightly lower than in the flat 
space case. Below we confirm this finding for NTVG conditions, with a very strong and 
complicated surface magnetic field. 

The gap electric field Ey along B s results from a deviation of the local charge density 
p ps from the corotational charge density p c = ((/a)pcj, where pcj = —ft • B s /(27rc) is 
the flat space-time Goldreich- Julian (1969) charge density, Q = 2-k/P, ( = 1 — n g , k 9 ~ 
(r g /R)(I /MR 2 ), a = (1 — r g j ' R) 1 ^ 2 is the redshift factor, r g is the gravitational radius, M 
is the neutron star mass and I is the neutron star moment of inertia. The potential V and 
electric fields E\\ within a gap are determined by GR-analog (Muslimov & Tsygan 1992) of 
the one dimensional the Poisson equation {1/ a)d 2 V / dz 2 = — 4ir(p — p c ) = 2QB s (/(ca), with 
a boundary condition E\\(z = h) = 0, where h is the height of an infinitesimal gap. The 
solution of the Poisson equation gives 

E\ l (z) = C(2n^.)(h-z), (1) 

and 

AV = (^h 2 . (2) 

In further calculations we will adopt a typical value of the correction factor ( ~ 0.85 (corre- 
sponding to M — 1M Q , R = 10 6 cm and / = 10 45 g cm 2 ), although its value can be as low 
as 0.73 (Zhang et al. 2000). Thus, the potential drop within the actual gap can be 15 to 27 
percent lower than in the conventional RS75 model. 
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The polar cap surface of the pulsar with ft ■ B < is heated by a back-flow of rel- 
ativists electrons (accelerated in the parallel electric field E») to the temperature T s = 
k}/ A {eAVN / aTrr 2 ) 1 / 4 , where AU is described by equation (2), N — ixr 2 p B s j eP is the Goldreich- 
Julian kinematic flux and the heat flow coefficient 0. 1 < A: < 1 is described in Appendix B. 
The thermal condition for the vacuum gap formation can be written in the form Tj/T s > 1, 
where T s is the actual surface temperature described above, and Tj = Ae c /30k is the iron 
critical temperature above which ^Fe ions are not bound on the surface (that is, a copious 
amounts of ^Fe ions at about Goldreich- Julian density are available for thermionic ejection 
from the surface), where Ae c is the cohesive energy of condensed ^Fe matter in the neu- 
tron star surface and k= 1.38 x 10~ 23 JK~ 1 is the Boltzman constant (Cheng & Ruderman 
1980; Usov & Melrose 1995). The properties of condensed matter in very strong magnetic 
fields characteristic of a neutron star surface have been investigated by many authors using 
different examination methods (for review see Usov & Melrose 1995). There exist many 
discrepancies in determination of the cohesive energy Ae c and in this paper we refer to 
the two papers, representing the limiting extreme cases. Abrahams & Shapiro (1991; AS91 
henceforth) estimated Ae c = 0.91 keV, 2.9 keV and 4.9 keV for B s = 10 12 G, 5 x 10 13 G 
and 10 13 G, respectively. These values were approximated by Usov & Melrose (1995) in the 
form Ae c ~ (0.9keV)(P s /10 12 G) a73 , which leads to critical temperatures 

Ti = (3.5 x 10 5 K)(P s /10 12 G) a73 = (6 x 10 5 )& a73 (PP_i 5 ) a36 K. (3) 

On the other hand, Jones (1986; J86 henceforth) obtained much lower cohesive energies 
Ae c = 0.29 keV, 0.60 keV and 0.92 keV for B s = 2 x 10 12 G, 5 x 10 12 G and 10 13 G, 
respectively. They can be approximated by Ae c ~ (0.18 keV)(P s /10 12 G) a7 and converted 
to critical temperatures 

Ti = (0.7 x 10 5 K)(P s /10 12 G) a7 = (1.2 x 10 5 )6 a7 (P • P_i 5 ) a36 K. (4) 

Below we consider the condition Ti/T s > 1, using both expressions for critical temperatures 
described by equations (3) and (4), for CR- and ICS-dominated NTVG models, separately. 

2.1. CR-NTVG 

In this case the gap height h = hen is determined by the condition that h = l p h, where 
l ph Rj sin 1Z = (B±/B S )TZ is the mean free path for pair production by a photon propagating 
at an angle 9 to the local surface magnetic field (RS75). The CR-NTVG model is described 
by the following parameters: the height of a quasi steady gap 

h CR = (3 x 10 3 )C- 3/7 ^ /7 6- 3 / 7 P 3 / 14 Pr i 3 / 14 cm, (5) 
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(notice typographical errors in eq. [6] of GM01, which are corrected here), the gap potential 
drop 

AV CR = (1.2 x lO^C^nfb^p-^P'Sl 4 V, (6) 
and the surface temperature 

T s = (3.4 x lO^C^k^Kfb^p-^Plil K. (7) 

The thermal condition Tj/T s > 1 for the formation of CR-NTVG leads to a family of critical 
lines on the P — P diagram (see Fig. 1 in GM01) 

P-is > A\ 016 k 114 K° 6 64 b- 2 p- 2 - 3 , (8) 

where A = (2.7 x 10 3 ) 1 / 2 = 52 for AS91 case (eq. [3]) and A = (3.96 x 10 6 ) 1 / 2 = 1990 for 
J86 case (eq. [4]). Alternatively, one can find a minimum required surface magnetic field 
B s = bB d expressed by the coefficient b in the form 

b c J n = AC om k^nr 2 p- 115 P:^. (9) 



2.2. ICS-NTVG 

In this case the gap height h = hies is determined by the condition h = l p h ~ l e , where l e 
is the mean free path of the electron to emit a photon with energy hw = 2 / ~fheB s /mc (Zhang 
et al. 2000, GM01). The ICS-NTVG model is described by the following parameters: the 
height of a quasi steady gap 

h ICS = (5 x lO 3 )^- 14 ^- - 07 ^ 57 ^ 1 ^" - 36 ^ 5 cm, (10) 
the gap potential drop 

AV ICS = (5.2 x lO 12 )^- 7 ^" - 14 ^- 14 ^ 1 ^" 1 - 22 ^ 5 V, (11) 
and the surface temperature 

T s = (4 x 10 6 )C ai8 A; a21 < 28 p- a43 K. (12) 

The thermal condition Tj/T s > 1 for the formation of ICS-NTVG leads to a family of critical 
lines on the P — P diagram (see Fig. 1 in GM01) 

P-i5 > B 2 c°- 5 k°- 7 n° 6 8 b- 2 p- 2 - 2 , (13) 
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where B = (2 x 10 2 ) 1 / 2 = 14 for AS91 case (eq. [3]) and B = (1.69 x 10 4 ) 1 / 2 = 130 for J86 case 
(eq. [4]). Alternatively, one can find a minimum required surface magnetic field B s = bB d 
expressed by the coefficient b in the form 

e = BC^'i^R^r ° 5 5 . (14) 

A comparison of spark developing time scales in ICS- and CR-dominated vacuum gaps is 
presented in Appendix C. 

2.3. NTVG development in pulsars with drifting subpulses 

GM01 examined the P — P diagram (their Fig. 1) with 538 pulsars from the Pulsar 
Catalog (Taylor et al. 1993) with respect to the possibility of VG formation by means of 
the condition Tj/T s , with T { determined according to AS91 (eq. [3]). They concluded that 
formation of VG is in principle possible, although it requires a very strong and curved surface 
magnetic field B s = b ■ B>d> 10 13 G, irrespective of the value of B d . Here we reexamine this 
problem by means of equations (9) and (14), with critical temperatures corresponding to 
both AS91 case (A = 52 and B = 14) and J86 case (A = 1990 and B = 130), for CR 
and ICS seed photons, respectively. Therefore, we cover practically almost the whole range 
of critical temperatures between the two limiting cases, which differ by a factor of five 
between each other. We adopt the GR-IFD correction factor ( = 0.85 (Harding & Muslimov 
1998; Zhang et al. 2000), the normalized radius of curvature of surface magnetic field lines 
0.01 < 1Z 6 < 1.0 (GM01 and references therein) and the heat flow coefficient 0.1 < k < 1.0 
(Appendix B). 

The results of calculations of NTVG models for 42 pulsars with drifting subpulses and/or 
periodic intensity modulations (after Rankin 1986) are presented in Fig. 1. We calculated 
the ratio B s /B q = 0.0453 6(PP_ 15 ) a5 , where B q = 4.414 x 10 13 G and the coefficient b is 
determined by equation (9) and equation (14) for CR and ICS seed photons, respectively. 
Noting that the functional dependence on P and P_is in both equations is almost identical, 
we sorted the pulsar names (shown on the horizontal axes) according to the increasing 
value of the ratio B s /B q (shown on the vertical axes). Calculations were carried out for 
B s /B q > 0.1, since below this value the near-threshold treatment is no longer relevant. On 
the other hand, the values of B s /B q are limited by the photon splitting threshold, which is 
roughly about 10 14 G in typical pulsars (Barring & Harding 1998, 2001; Zhang 2001, see 
also astro-ph/0 102097). Therefore, the physically plausible NTVG models lie within the 
shaded areas, with the upper boundary determined by the photon-splitting level and the 
lower boundary determined by the lowest ICS line calculated for k — 0.1 (see below). The 
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left hand side of Fig. 1 corresponds to AS91 case (A = 52 and B = 14 in eqs. [9] and 
[14], respectively) and the right hand side of Fig. 1 corresponds to J86 case (A = 1930 and 
B = 130 in eqs. [9] and [14], respectively). Four panels in each side of the figure correspond 
to different values of the radius of curvature TZq = 1.0, 0.1. 0.05 and 0.01 from top to bottom, 
respectively (indicated in the upper corner of each panel). Two sets of curved lines in each 
panel correspond to the CR (thin upper lines) and ICS (thick lower lines) seed photons, 
respectively. Three different lines within each set correspond to different values of the heat 
flow coefficient k = 1.0 (dotted), k = 0.6 (dashed) and k = 0.2 (long dashed), and lower 
boundary of shaded areas corresponds to k — 0.1. 

A visual inspection of the model curves within shaded areas in Fig. 1 shows that in AS91 
case ICS-NTVG is favored for larger radii of curvature 1Z 6 > 0.05, while CR-NTVG requires 
lower values of 1Z e < 0.1. In J86 case, only ICS-NTVG corresponding to TZ& < 0.1 can 
develop. If the actual cohesive energies correspond to some intermediate case between AS91 
and J86 cases, they will also be associated with ICS-NTVG. Therefore, we can generally 
conclude that ICS-NTVG model is apparently favored in pulsars with drifting subpulses. 
CR-NTVG is also possible, although it requires a relatively low cohesive energies, as well as 
an extremely strong [B s jB q ~ 1) and/or curved (7Z 6 ~ 0.01) surface magnetic field at the 
polar cap. 



3. Conclusions 

There is a growing evidence that the radio emission of pulsars with systematically drift- 
ing subpulses (grazing cuts of the line-of-sight) or periodic intensity modulations (central 
cuts of the line-of-sight) is based on the inner vacuum gap developed just above the polar 
cap (Deshpande & Rankin 1999, 2001; Vivekanand & Joshi 1999; Gil & Sendyk 2000). To 
overcome the binding energy problem Xu et al. (1999, 2001) put forward an attractive but 
exotic conjecture that pulsars showing the drifting subpulses represent bare polar cap strange 
stars (BPCSS) rather than neutron stars. However, as demonstrated in this paper, invoking 
the BPCSS conjecture is not necessary to explain the drifting subpulse phenomenon. The 
quasi steady vacuum gap, with either curvature radiation or inverse Compton scattering seed 
photons, can form in pulsars with O ■ B < 0, provided that the actual surface magnetic field 
at the polar cap is extremely strong B s ~ 10 13 G and curved 1Z < 10 6 cm, irrespective of the 
value of dipolar component measured from the pulsar spindown rate. We have used two sets 
of the cohesive (bounding) energies of the surface iron ions: higher values obtained by AS91 
and about five times lower values obtained by J86. If the actual cohesive energies are close 
to J86 values, then only ICS controlled VG can form and CR controlled VG never forms 
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even under the most extreme conditions (see Fig. 1). 

It is worth noting that ICS discussed in this paper is the so-called "resonant" scattering. 
The choice of this mode of ICS is justified by the superstrong magnetic fields relevant for 
the near-threshold regime. However, it is clear that the thermal photons also contribute 
some pair-producing gamma-rays (e.g. Zhang et al. 1997). This process cannot dominate 
the breakdown of the gap, since it requires higher surface temperatures, which may inhibit 
the formation of VG in the first place. One can only speculate that the dominating thermal 
ICS mode is associated with the pulse nulling phenomenon. These issues require further 
investigation. 

The pulsars with drifting subpulses and/or periodic intensity modulation do not seem 
to occupy any particular region of the P — P diagram (see Fig. 1 in GM01). Rather, 
they are spread uniformly all over the P — P space, at least for typical pulsars (excluding 
young and millisecond pulsars, in which observations of single pulses are difficult in the first 
place). Therefore, it seems tempting to propose that radio emission of all pulsars should 
be driven by vacuum gap activities. An attractive property of such proposition is that 
the nonstationary sparking discharges induce the two-stream instabilities that develop at 
relatively low altitudes (Usov 1987; Asseo & Melikidze 1998) where the pulsars radio emission 
is expected to originate (Cordes 1978; Kijak & Gil 1997, 1998). It is generally believed that 
the high frequency plasma waves generated by the two-stream instabilities can be converted 
into coherent electromagnetic radiation at pulsar radio wavelengths (e.g. Melikidze et al. 
2000). With such scenario, all radio pulsars would require a strong, non-dipolar surface 
magnetic field at their polar caps (e.g. Gil et al. 2002 a, b). 

This paper is supported in part by the KBN Grant 2 P03D 008 19 of the Polish State 
Committee for Scientific Research. We are indebted to anonymous referee for insightful 
comments and constructive criticism, which greatly helped to improve the final version of 
the paper. We are grateful to Prof. V. Usov and Dr. B. Zhang for helpful discussions. We 
also thank E. Gil, U. Maciejewska and Dr. M. Sendyk for technical help. 

A. Field emission 

The vacuum gap can form in pulsars with Q ■ B < if the actual surface tempera- 
ture T s is not high enough to liberate ^Fe ions from the polar cap surface by means of 
thermal emission. Now we examine the field (cold cathode) emission, which is possibly 
important when thermionic emission is negligible (e.g. Zhang & Harding 2000). The max- 
imum electric field (along B s ) at the NS surface E\\(max) = ((4-K/cP)B s h IC s = (1-25 x 
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10 9)£0.86 6 -i 7 jo.57_ p o.i4 ^ / cm, or taking into account that C = 0.85 and b>2{P ■ P_i 5 )" - 5 , 
E\\(max) < (5 x 10 8 )^- 57 P a64 P° 1 5 5 V/cm. The critical electric field needed to pull ^Fe ions 
from the NS surface is EAcrit) = (8 x 10 12 )(Ae c /26 keV) 3 / 2 V/cm (Abrahams & Shapiro 
1991; Usov & Melrose 1995, 1996), where the cohesive (binding) energy of iron ions in a 
strong surface magnetic field B s ~ 10 13 G is Ae c = 4.85 keV (Abrahams & Shapiro 1991). 
Thus, E\\(crit) = 6.4 x 10 11 V/cm >> E\\ ~ 10 9 V/cm and no field emission occurs. It is 
possible that the cohesive energy is largely overestimated and Ae c can be much smaller than 
about 4 keV even at B s > 10 13 G. We can thus ask about the minimum value of Ae c at 
which the field emission is still negligible. By direct comparison of E\\(max) and E\\{crit) 
we obtain Ae c > 40a: - 67 eV, where x = 7^- 57 P a64 P°f 5 is of the order of unity. Therefore, 
contrary to the conclusion of Jessner et al. (2001), no field emission is expected under any 
circumstances. 



B. Heat flow conditions at the polar cap surface 

Let us consider the heat flow conditions within the uppermost surface layer of the pulsar 
polar cap above which NTVG can develop. The basic heat flow equation is (e.g. Eichler & 
Cheng 1989; Usov 2001) 

C^L = A ( K f*E\ (Bl) 
dt dx \ " dx J ' 

where C is the heat capacity (per unit volume) and km >> k± ~ is the thermal conductivity 
along (||) surface magnetic field lines (which are assumed to be perpendicular (_L) to the polar 
cap surface). The heating of the surface layer of thickness Ax is due to sparking avalanche 
with a characteristic development time scale At. We can write approximately dT/dt xs 
T/At, dT/dx « T/Ax, d/dx(dT/dx) w T/Ax 2 and thus C/At w k\\Ax 2 . Therefore, 
the crust thickness Ax heated during At is approximately Ax ss (/t||At/C) (within an 
uncertainty factor 5x of 2 or so), where the time scale corresponds to the spark development 
time At ps 10/is (RS75). The energy balance equation is Qheat = Qrad + Qcond, where 
Qheat = encjAVmax (e.g. RS75), Qrad = (T s is the actual surface temperature and 
a = 5.67 x 10~ 5 erg cm _2 i^ _4 s _1 ), and Q con d — —K,\\dT/dx ~ — KnT s /Ax. We can now 
define the heat flow coefficient k = Q ra d/ Qheat < 1, which describes deviations from the 
black-body conditions on the surface of the sparking polar cap. In other words, the value of 
k describes the amount of heat conducted beneath the polar cap which cannot be transferred 
back to the surface from the penetration depth Ax during the time-scale At (see Appendix 
C). The coefficient k can be written in the form 

u — \ — \ (no) 

1 + kJ(oT*Ax) l + {K ll C) 1 / 2 /(aT^At 1 / 2 y { ' 
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whose value can be estimated once the parameters C, k\\ and T s as well as Ax or At are 
known. 

The matter density at the neutron star surface composed mainly of ^Fe ions (e.g. Usov 
& Melrose 1995) is p(B s ) - 4x 10 3 (5 S /10 12 G) 6 / 5 g cm" 3 (Flowers et al. 1977). Thus for 
B s ~ (1 3) x 10 13 G we have p ~ (0.6 -i- 2.4) x 10 5 g cm -3 . The thermal energy density of 
^Fe is U Fe ~ 2.2 x lO^Tg erg cm" 3 (Eichler & Cheng 1989), where p 5 = p/(10 5 g cm" 3 ) 
and T 8 = T s /(10 8 K). Thus, the heat capacity C = dU Fe /dt = lO~ 8 dE Fe /dT 8 = (1 + 5) x 
10 n p5 erg cm" 3 K _1 . For p 5 « 1 we have C ~ 2x 10 11 erg cm _3 K" 1 . The longitudinal thermal 
conductivity can be estimated as /ty ~ (2^-4) x 10 13 erg s _1 cm _1 K _1 (see Fig. 5 in Potekhin 
1999). Thus K\\/C ~ 1.5 x 10 2 cm 2 and the penetration depth Ax p=j 10At 1 / 2 s" 1 / 2 cm. Since 
the characteristic spark development time scale At ~ lOps = 10~ 5 s (see Appendix C), then 
Ax 0.03 cm. More generally, the penetration depth can be written as 



where k 13 = k\\/10 13 2-h4 and Cn = C/10 11 ^ 1^-5 and the uncertainty factor 5x ps 0.5^2. 
Now the heat flow coefficient can be written as 



where T 6 = T s /10 6 and A = (n n Cu) 1/2 / 5x ps 0.7 + 6.3. Figure 2 shows variations of the 
heat flow coefficient k versus the surface temperature T 6 (in 10 6 K) for three values of A = 1 
(upper curve), A = 3 (middle curve) and A = 6 (lower curve). As one can see, for realistic 
surface temperatures of a few times 10 6 K, the values of the heat flow coefficient k are in 
the range 0.2 0.9. Gil et al. (2002a) found from independent considerations that in PSR 
B0943+10 the heat flow coefficient k < 0.8, in consistency with the results obtained in this 
paper. 

The above treatment of the heat transportation corresponds to one spark event or a 
relatively short sequence of sparks reappearing at the same place of the polar cap (modulo 
the slow E x B drift). If a longer sequence of consecutive sparks can operate at the same 
place, deep-layer photons can also diffuse out of the surface, although after a considerably 
longer time (Eichler & Cheng 1989). The contribution from this component would make k 
closer to unity. Even if this is the case, ICS driven VG is still possible to form for all pulsars 
from our sample, although an extremely strong and curved surface magnetic field would be 
required (the right hand-side panels of Fig. 1). 




(B3) 
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C. Characteristic spark development time scale 



The characteristic spark development time scale is defined as the time interval af- 
ter which the density of electron-positron spark plasma grows from essentially zero to 
the Goldreich- Julian (1969) value. In the case of CR-induced sparks this time scale is 
t cr — (30 — A0)h C R/c ~ 10/is. The time structure of the ICS-induced spark avalanche 
may be different. This follows from the fact that in the CR-induced avalanche the mean 
free path l e of electron/positron to generate one pair-producing gamma photon is much 
smaller than the mean free path l p h of gamma photon to produce an electron-positron pair, 
while in the ICS-induced avalanche l e is comparable with l ph (Zhang et al. 2000). Let us 
introduce two time scales: the time t p = l e /c for a charged particle (electron or positron) 
to emit a gamma- quantum, and the time t q = l p h/c for a gamma-quantum to create an 
electron-positron pair. In both CR- and ICS-induced gap discharges the condition t p < t q 
holds (t p <C t q in the CR case and t p ~ t q in the ICS case), which suggests that the 
avalanche time scale is determined by t q . In fact, one can assume that the cascade con- 
sists of a sequence of n(t) ~ t/t q generations, occurring between the initial time t — 
and the current time t. Thus, the number of particles created during the n-th genera- 
tion is N(t) = (2t q /t p ) n ^ = (2t q /t p ) t f tq . This estimate represents the following scheme of 
generation-by-generation cascade development: (2t q /t p ) — > (2t q /t p ) 2 — > ... — > (2t q /t p ) t ^ tq . 
For the ICS-induced cascade we have t p ~ t q ~ hjcs/c, while for the CR-induced case 
tq = hcn/c and t p ~ 2^v m hj 'Lqr ~ 3 x 10~ 3 7?.6/7, where Lqr = (2/3)ce 2 7 4 /7?. is the power 
of the curvature radiation at the characteristic frequency v m = OA^c/R k 2 x 10 9 7 3 /7?.s~ 1 , 
7 = eAV/mc 2 = 10 7 C°' 14 7?.g 57 6 0,14 is the Lorentz factor of relativistic charged particles, 
AV = AV C r is the NTVG potential drop described by eq.(6), and 1Z = 7£ 6 10 6 cm is the 
radius of curvature of the surface magnetic field B s = bB d (we assumed P = P_i 5 = 1 for 
simplicity). Here e is the electron charge, c is the speed of light and h is the Planck constant. 

The total number of particles created during the time interval t is iV = J* N(x)dx, 
where x = t/t q . Thus, for the CR-induced case in which t q /t p ~ Hcr/ (t p c) 




(CI) 



and for the ICS-induced case in which t q /t p ss 1 



tc 



Nws 




(C2) 



In 2 



where hcR and hies are described in eqs.(5) and (10), respectively. 
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The cascade terminates at the time t — r when the spark plasma density reaches the 
Goldreich- Julian density n GJl that is N IC s/h 3 ICS = N C R/h CR ~ n GJ . Using this condition 
and eqs.(Cl) and (C2) we obtain the ICS-induced spark development time scale tjcs as a 
function of CR-induced spark development time scale tqr in the form 

hies , f h cn\ . hies ■> I ln2 fh IC s\ 3 \ 

Tlcs= >^ u { 2 i^) Tc « + ^(^^{h^) )■ (C3) 

This time scale depends also on other pulsar parameters, mostly on b and TZq (and weakly 
on P,P,k and (). The second term in eq. (C3) is always negative and much smaller than 
the first term (for t C r > 201icr/c; e.g. RS75). Therefore we obtain an upper limit for the 
ratio 

^. 2 . 4c o, V „ W , 9 ^_«^). (C4) 

Fig. 3 shows the upper limits (equality in eq. [C4]) of the ratio tics/ t cr as a function of 
B s /B q = 0.04536 (see section 2.3 and Fig. 1 for explanation), for different values of 1Z 6 
and k (for the clarity of presentation we adopted P = P_i 5 = 1 and ( = 0.85). As one 
can see from this figure, the characteristic spark development time scales of CR- and ICS- 
induced discharges are comparable. This is a consequence of the fact that the stronger the 
magnetic field the smaller the l e for ICS, thus the larger the pair-production multiplicity 
and the smaller the spark development time scale Tics- We can thus generally state that in 
the near-threshold regime ICS-driven spark avalanche looks very similar to the conventional 
RS75 CR-driven avalanche. 
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Figure Captions 



Fig. 1. — Models of NTVG for 42 pulsar with drifting subpulses and/or periodic intensity- 
modulations for AS91 and J86 cohesive energy calculations. Both curvature radiation (CR) 
and resonant inverse Compton radiation (ICS) seed photons generating electron-positron 
pairs in a superstrong surface magnetic field B s are considered. The ratio of B s /B q > 
0.1 (near threshold treatment ), where B s is the actual surface magnetic field and B q = 
4.4 x 10 13 G is the quantum critical magnetic field, is presented on the vertical axis. The 
horizontal axis corresponds to 42 pulsars with drifting subpulses and/or periodic intensity 
modulations. Three different lines correspond to different values of the heat flow coefficient 
k = 1.0 (dotted), k = 0.6 (dashed) and k = 0.2) (long dashed). Four different panels 
correspond to different values of curvature radii from 1Z e = 1.0 (top) to 1Z 6 = 0.01 (bottom). 
The physically plausible models correspond to shaded areas limited by the photon splitting 
threshold B s ~ 10 14 G from above and by the lowest ICS curve (calculated for k = 0.1) 
from below. 

Fig. 2. — Heat flow coefficient k versus surface temperature T 6 = T s /10 6 K (see text for 
explanations). 

Fig. 3. — Upper limits for the ratio of development time scales in ICS- and CR- induced 
spark discharges of NTVG, as a function of the surface magnetic field B s /B g , ranging from 
0.1 (near threshold treatment) to 2.25 (photon splitting limit). Four values of TZ e =1.0, 0.1, 
0.05 and 0.01 (from the top to the bottom) and two values of k—1.0 (solid line) and 0.2 
(dashed line) were used. 
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